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Abstract 

This work proposes an analytical framework to study how relay selection strategies perform in 
half- and full-duplex deployments by combining renewal theory and stochastic geometry. Specifically, 
we assume that the network nodes - operating in either half- or full-duplex mode - are scattered 
according to a two-dimensional homogeneous Poisson point process to compute the relay selection 
cost by using a semi-Markov process. Our results show: (i) fixed relay outperforms the reactive option 
in either cases, (ii) the performance of both reactive and fixed relay strategies depends on the self¬ 
interference attenuation in full-duplex scenarios, evincing when they outperform the half-duplex option, 
and (in) the reactive relay selection suffers from selecting relays at hop basis, while the fixed relay 
selection benefits most from the full-duplex communication. 


I. Introduction 

In the recent years, Full-Duplex (FD) communication has gained considerable attention from 
both academy and industry (T), [[2]|. By allowing simultaneous transmission and reception on the 
same frequency band, FD networks can potentially double the spectral efficiency compared to 
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current Half-Duplex (HD) schemes. However, such potential is harmed by Self-Interference (SI) 
from the transmit to receive antenna 03, 03 • Due to recent advances in antenna design associated 
with analog and digital interference cancellation, most of the SI can be mitigated, and therefore 
FD communication becomes feasible and are a step towards meeting high demands of spectral 
efficiency of 5G 0. 

Current works have shown that FD is a viable solution for small cell deployments [0, due 
to not only cost and size constraints but also to compatibility with current HD systems. Then, 
due to the FD capability, the base station can, for instance, simultaneously schedule uplink and 
downlink transmissions, see for instance 0. Moreover, base station can also act as relays for 
the legacy network as well as for its own subscribed users. Cooperative diversity appeared as a 
way to combat the effects of the fading, by allowing a single antenna user to experience spatial 
diversity 0. Such schemes perform even better when the relay operates in a FD fashion 0. 

Bearing this in mind, and in the context of a dense deployment of small cells, relaying 
selection become an attractive solution given a trade-off between complexity and efficiency 
0, Q. Relay selection was initially proposed in Q for HD cooperative networks and has 
gained considerable attention since. Nonetheless, most of the works so far have focused on 
HD selection algorithms, except for 0, 0 where the authors focus on the performance of a 
three-node relaying scenario. For instance, 0 focus on selection algorithms for FD amplify- 
and-forward and model the residual SI at relay as a Rayleigh random variable. While in 0, 
the authors address the performance of a FD decode-and-forward protocol and assume simpler 
residual SI model (constant attenuation factor). 

In this work, we assess the performance at the network level of a relay selection procedure 
for cooperative diversity protocols in FD mode. To do so, we model the relaying selection 
procedure as a semi-Markov process and then we investigate the impact of the SI and Co- 
Channel Interference (CCI) on the network performance. Moreover, a dynamic relay selection 
procedure is considered where a suitable relay is selected at each hop and the cost of this selection 
procedure is incorporated into the achievable rate. And yet, we resort to stochastic geometry to 
model network deployments and capture the dynamics of the network, more specifically the 
aggregate interference at the receiver, which can be characterized through a cumulant-based 
framework introduced in m, cn. 

Our main contributions are summarized as follows: 


3 


• The standard binary tree algorithm is used to implement a relay selection procedure for FD 
networks; 

• The framework originally introduced in ||T2| is extended to account for the FD relaying 
operation; 

• The performance of the aforesaid relay selection protocol are evaluated and compared to 
its HD counterpart; 

• Different from |8|, our network level analysis accounts for a more detailed model of the 
residual SI based on lfl3l and also characterizes the CCI at the receiver. Actually, instead of 
modeling the residual SI by a constant value, we consider a more detailed model through 
a random variate drew from a Ricean distribution with large K -factor. 

The remainder of this paper is organized as follows: Section HH describes the problem under 
investigation and introduces the system model. Section [HD describes the network topology 
and introduces our analytical framework based on the stochastic geometry, while Section HVl 
characterizes the relay selection procedures and the probability generating function concept. 
Section |V] introduces the semi-Markov process that is used to evaluate the performance of the 
proposed relay selection mechanism. In Section ED performance metrics are introduced and then 
utilized to obtain the numerical results and carry out their discussion. We draw conclusions and 
final remarks in Section IVII1 

II. Problem Description and System model 

As aforesaid, full duplex networks constitute a promising alternative to achieve high data 
rates in the upcoming 5G systems. In such networks, cooperative nodes benefit most from the 
inherent rate increase allowed by transmitting and receiving simultaneously. Typically, available 
literature concentrates on assessing the throughput that is attainable by combining full duplex 
and cooperative diversity protocols, and overlooks the impact of the underlying relay selection 
on the final achievable data rate J8J. 

In this paper, we contribute to previous results by introducing an analytical framework that 
combines stochastic geometry and semi-Markov processes to account for the cost of selecting a 
suitable relay. On that regard, we initially use stochastic geometry to derive the distributions of the 
received power and aggregate interference experienced by the receiver of interest in Section [Till 
Thereafter, we use the Probability Generating Function (PGF) concept lfl4l . lfI51 to characterize 
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the Contention Resolution Interval (CRI) (holding time) of the relay selection algorithms as can 
be seen in Sections [TV] and [Vj Our investigations extend lfl2l in which Authors introduce an 
analytic framework based on semi-Markov processes with constant holding time and apply it to 
study an one-dimensional network model. 

We assess the performance of interference networks operating in FD mode under composite 
fading channel - Lognormal shadowing and Nakagami-m fading. The Downlink (DL) of a 
traditional HD network constitutes our benchmark scenario wherein the user of interest is 
interfered by surrounding small cells. Base stations independently schedule a random user 
terminal in every transmission interval. All communicating nodes use antennas with omni¬ 
directional radiation pattern. Base stations and User Equipments (UEs) are also assumed to 
have full buffer and symmetric traffic patterns Ifl6ll . 

Consider all channels are quasi-static and undergo a composite fading distribution [fT71 . which 
is approximated by a lognormal distribution with mean and variance (in logarithmic scale) given 
by /idB = £ [ip (m ) — In (m)] + and ct 2 b = (2, m ) + cr^ p , where ip (m) is the Euler psi 

function and ((2. m) is the generalized Riemann zeta function |fl8ll . Radio links are affected 
by path-loss, large-scale shadowing and multi-path fading which are assumed to be mutually 
independent and multiplicative phenomena |[T9l . The received power at the user of interest u 0 
from an arbitrary transmitter bi located r ?:o meters away is Y i0 = p i0 , where p i0 yields the 

transmit power, a is the path-loss exponent and :r,o represents the composite fading channel. In 
what follows, when there is no fear of ambiguity let us omit the subscripts of the interference 
component terms. Fig. Q] illustrates a realization of the random network topology where UEs 
and small cell Base Stations (BSs) are uniformly scattered over network area. 

III. Network deployment and Stochastic Geometry 

Active nodes constitute a homogeneous Poisson Point Process (PPP) $ with intensity A in M 2 . 
The number of nodes in an arbitrary region 1Z of area A is a Poisson Random Variable (RV) 
with parameter A A KOI . 

The fading effect is assumed as a random mark associated with each point of $, thus 
<f> = {(p,x) ;ip E <£>} corresponds to a marked point process on the product space M 2 x M + , 
whose random points <p denoting transmitters locations and belong to $ f20l . [f2TTl . The network 
deployment model is given by a spatial PPPs $ BS (T UE ), whose random points <p represent 
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Fig. 1. Illustration of a random (bi-dimensional) deployment of base stations and user terminals over an arbitrary network area. 
Shaded squares represent user terminals, while shaded up-triangles depict base stations. A heat map represents the corresponding 
random composite fading channel where the fading intensity varies from red/strong to dark blue/weak. 


the locations of BSs (UEs). In addition, we consider a fixed number of potential relays that 
detect source’s message and take part in the selection procedure, it is well known that such 
conditional distribution of points within an arbitrary forwarding region A is described by a 
binomial point process lf2Til . As shown next, by employing stochastic geometry concepts to this 
deployment model, it is possible to conveniently derive closed form expressions to characterize 
the distributions of the received power and corresponding aggregate CCI with respect to the 
receiver of interest. Indeed, the path loss attenuation (function of the distance to the tagged 
user), the lognormal shadowing and Nakagami-m fading are combined using the framework of 
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urn m to obtain the distribution of the received power at the user of interest. 

A. Received Power from a Random Transmitter 

We resort to the framework introduced in ifTOl . fTTI to attain the characteristic function and the 
respective cumulants (k) of the received power and aggregate interference of random transmitters 
within the forwarding region of the user of interest. Consider a RV Y = R~ a X describing the 
power received at the tagged receiver from a random transmitter in O (as defined in Section UT]), 
with R varying from R m to R M , and X ~ Lognormal(// (1 B. <t<ib) representing the squared- 
envelope of the composite fading channel. Then, the characteristic function of Y is [11, Appendix 
B] 'ky (cc) = 2(R 2 m — Rffj^Ex [R (cu)] where R(cc) = J'^ M exp (jcopr~ a x)rdr and E x [ • ] 
yields the expectation of the enclosed expression over the RV X, and is a general formulation 
which characterizes the distribution of any random transmitter within the reception range of a 
tagged receiver. Thereafter, the nth cumulant of Y is given by [QTj, Appendix C] 

1 n 

Kn = ~ y(A ))-B n ,k [/?1, /?2, • • • , An-fc+l)] , (1) 

X 

k=0 

where i yields the imaginary unity, g ( u ) = In (n), B n>k [/3 1: /3 2l ..., is the partial Bell 

p2-na _ jj2—na 

polynomial flU and 8 n = i n p n x E x [x n \. 

na — 2 

Next, we use the aforesaid cumulant-based framework to characterize the aggregate inter¬ 
ference at the user of interest for both the HD and FD configuration scenarios. Under the 
assumptions of Section HO and with respect to the tagged receiver, we now use our analytical 
framework to characterize the resulting aggregate CCI in two Evaluation Scenarios (ESs): ESI 
is the benchmark scenario where BSs and UEs operate in HD mode; while in ES2 both BSs 
and UEs operate in FD mode. 

B. Aggregate CCI from a Poisson Field of Interferers 

In our benchmark scenario, denoted HD relaying, all nodes operate in HD mode. The corre¬ 
sponding aggregate CCI at a tagged receiver is then given by Zq D = Si )ei BS ^o- Notice 

that the tagged receiver only experiences interference coming from nearby base stations since 
we assume the DL as the benchmark scenario. 
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In this context the aggregate CCI at the tagged receiver is IflOl Section V] 

K ”( i ' BS ) = IrS Ell °- to] ■ <2) 

C. Aggregate CCI in the FD Configuration 

In the FD configuration, the tagged receiver is subject to the interference from FD serving base 
stations, other FD UEs and its intrinsic SI component. The self-interference channel between 
transmitting and receiving antennas of a full duplex transceiver exhibits a strong line of sight 
component and can be represented using a Ricean distribution with large A"-factor ||24ll (a K- 
factor of 14.8dB is used, which corresponds to a Nakagami distribution with the parameter 
m = 16). As a result, the aggregate CCI at the tagged receiver is characterized by 

Z f 0 D = Sp 00 x 00 + Y i0 + Yj 0 , (3) 

(<pj,Xj)e® UE 

where p 00 and 5 represent the SI component and the corresponding attenuation factor, respec¬ 
tively. 

Since BSs and UEs are assumed to be independently scattered over the network deployment 
area, the resulting processes (<I> BS and $ UE , respectively) from each such tier are also independent 
lf20l . Therefore, the cumulants additivity property is employed to obtain the aggregate CCI from 
multiple interfering tiers ifTOl . In other words, if U and V are independent RVs, then, we can 
derive the cumulant of their sum as n n (U + V) — n n (U) + K n (V). 

D. Signal-to-Interference Ratio and Outage Probability 

By using the results from previous sections, it is straightforward to compute the Signal-to- 
Interference Ratio (SIR) and outage probability of the evaluation scenarios under study. The SIR 
and outage probability are used to assess how the relay selection strategies perform in HD and 
FD configuration mode as follows. 

Proposition 1: Let V 0 and V be Normal RVs (in logarithmic scale) representing the power 
received from the desired transmitter and the CCI at the receiver of interest. Under the assumption 
of the composite fading with Gamma-lognormal distribution, the SIR at the tagged receiver 
is T ~ Normal (/iy 0 — p v , Oy Q + Oy), and the outage probability is given by Pr [r < 7 th ] = 
Q [(/i r — 7 th ) /exp] , where /ip = pv 0 ~Fv an d of = ay + cry. Notice that the lognormal parameters 



can be attained as follows /j = In yn\/yj n\ + k 2 j, and cr 2 = In (1 + k 2 /h\), where K\ and n 2 
are given by (HJ. 

Proof: The SIR distribution is given by the quotient of two independent lognormal RVs, 
namely, e Vo which is the received power from the target transmitter, and e v which is an equivalent 
lognormal RV approximating the aggregate CCI at the tagged receiver ll25l . ■ 

IV. Relay Selection Algorithm and Probability Generating Function 

Hereafter, we characterize the cost of selecting a suitable relay at hop basis on the achievable 
throughput. To do that, we first account for the relay selection overhead using the PGF of the CRI 
and later combine the impact of the network characteristics with the relay selection mechanism 
throughout a semi-Markov process. In this work, a totally random approach based solely on 
the standard tree algorithm is used to implement the relay selection algorithm Il26l . It is worth 
noticing we use the conditional PGF of the CRI length ( L N ) that is derived in lfl4l . llT5ll for 
the same standard tree algorithm under investigation. Next, due to space limitations, we only 
summarize the main results regarding the relay selection procedure. 

A. Relay Selection Algorithm 

Two relay selection approaches are considered in this work: (i) a static one whereby a relay 
node is preassigned to forward source packets; and (ii) a reactive strategy by which a relay node 
is identified based on the selection criteria - i.e., the neighboring node that provides the longest 
advancement towards the final destination at each hop. While the former configuration may be 
more applicable to legacy networks where a smallcell may provide the last hop to cell edge 
users, the latter is tailored to upcoming 5G networks where humans and machines communicate 
in large scale deployment scenarios ll27l . On one hand, the static relay procedure relies on long 
term characteristics of the deployment scenarios, i.e. the relay is selected a priori following a 
network planning (it is arbitrary in our investigations) and does not change. On the other hand, 
the reactive relay procedure selects a random (new) relay within the transmission range every 
time a packet needs to be forwarded to the destination [fTIil . [fl4l . Il28ll . 

Succinctly, the reactive relay selection algorithm is implemented through a carrier sense 
collision avoidance mechanism similar to the IEEE 802.11 handshake using Request to Send 
(RTS)/Clear To Send (CTS) messages. The source transmits a RTS packet to initiate the relay 



selection procedure. Nodes that listen to this request reply with a CTS packet based on the 
predetermined probability of accessing the channel, p. If a collision occurs, nodes that have 
transmitted in previous slot retransmit or not based on a random access process similar to a Q- 
sided coin. To select the best suitable relay (greed criteria), i.e. the closest node to the destination 
whether there is one available, the source node should receive the replies from all the candidate 
relays [|T4|. Ifl5l . We consider the conditional CRI length when N nodes initially collide (which 
corresponds to a binomial point process ETl ). Here, aiming at shorten the selection procedure^, 
the first node (at a random location within transmission range) to successfully reply to the source 
is selected as the next hop relay - note that this mechanism differs from the one used in |fl4|. 
Ifl5l which waits for the replies of all neighboring relay candidates in order to select the most 
suitable based on the predefined selection criteria. 

We use the relay selection mechanism introduced in |[T4|. llT5l which combines the standard 
tree algorithm to resolve contentions and sectoral decision regions to identify the potential relays. 
Fig. [2] illustrates a snapshot of the sectoral decision region that is used in conjunction with the 
standard tree-based solution. In the next section, the PGF is used to characterize the distribution 
of the CRI associated with the selection of a random relay at hop basis. 

Remark 1: Notice that current 802.11 standards implements the RTS/CTS handshake as an 
optional feature to avoid frames collisions; besides, more recent standards such as 802.22 and 
802.11af also incorporate cognitive radio functionalities which makes it even easier to deal with 
side-information Ifl4l . 

Remark 2: In LTE systems, the contention-based random channel is based on ALOHA- 
type communication and is susceptible to performance degradation due to a high collision 
probability in the transmission of the preambles in dense deployments. Indeed, the contention- 
based operation of the LTE Physical Random Access Channel is also established by the signaling 
exchange of messages, namely, preamble transmission, random access response, connection 
request and contention resolution If27l . 11291 . 


’At the cost of selecting a relay that is not necessarily the closest one to the final destination 
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dest. 


Fig. 2. Illustration of the sectoral decision region with Q = 2 splitting groups and angular aperture of 180 degrees. Dashed 
lines defined the forwarding region (shaded regions in light gray and blue), black circles identify awake nodes and white circles 
identify asleep nodes. All awake neighbors within the shaded region are eligible relays. 


B. PGF of the Relay Selection Algorithm 

The conditional CRI length considering a Q-sidcd fair coin is generalized as, 

[ 1 , if N = 0,1; 

Ljm — < Q (4) 

1+ E Li,, if N >2, 

l 3=1 

where Ij is a discrete RV describing the number of candidate relays that tossed the j value of 
the Q-sidc coin. For a binary splitting tree, the PGF of the time spent by a random node in the 
relay selection procedure is give by the following proposition. 

Proposition 2: Under the assumptions of a binary standard tree algorithm (unbiased coin) and 


































11 


N conflicting nodes, the PGF of the CRI for a tagged packet is 
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(z), 

(7) 


k =0 


where Gi = z, G $ +1 (z) accounts for the CRI when the tagged node joins the first subset 
(flipped 0), x (z) describes the resolution interval when the tagged node joins the second 
subset (flipped 1), and B Nk = Pr [k flipped 0 | N contending nodes]. Note that the tagged packet 
may be in each subset with equal probability. 

Proof: Let N > 1 be the number of contending nodes that start the CRI. For a fair binary 
splinting tree we have B Nk = ( J ^)2 -fc . Thus, if the first splitting group has k members, the 
second contains the remaining N — k nodes. Each such sub-CRI is statically indistinguishable 
from a general CRI with equal number of elements. Therefore, the recursive formula Qn(z) = 
Yl!k=oB N ,kQk{z)Q n -k{z) is obtained. By the same line of reasoning, we can defined 

Gn+i{z) = Pr[f 0 = k\N other nodes collide] z k ; (8) 

G$ + i(z) = Pr[f 0 = k\N other nodes collide 

and the tagged node flips s]z k , (9) 

where to yields the time a tagged packet spends in the CRI before being successfully transmitted. 
Finally, computing the corresponding probabilities for a fair binary spiting tree, we obtain ©. 

■ 

In other words, Proposition [2] gives the time a tagged packet takes to be successfully transmit¬ 
ted. When a random node within source’s transmission range is selected, its distance distribution 
is characterized by CO)- To determine the holding time matrix, we resort to the factorial moment 
concept so as to compute the mean and variance of CRI for a tagged packet from Proposition 
©. 

Definition 1: Let X be a discrete RV taking non-negative integer values; then, the k\h factorial 
moment of X is given as E [X(X — 1) • • • (X — k + 1)] = GP0(1“), k > 0 ll30l . 
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Then, we obtain the following results from Definition Q] as G(l) = Pr [X = k], G'( 1) = E[X] 

and G"( 1) = E[X(X - 1)]. 


V. Semi-Markov Process 


Markov processes permit represent random systems whose outcome at any given instant 
depends only on the outcome that proceeds it (memoryless property, if the present is specified, the 
past has no influence in the future). A Markov chain is special kind of Markov process wherein 
the future evolution of the process depends on the present state and not on how it arrived at that 
state Oil . Typically, such Markov models have the property that a transition occurs at every 
instant - even if the transition returns to the previous occupied state, it occurs anyway. In this 
work, we use a semi-Markov process that is associated with the embedded Markov chain in 
Fig. [3] that tracks the relay selection protocols states. The semi-Markov process theory allow us 
to obtain the steady-state average delay and throughput during a cycle taking any state of the 
chain as a reference ll30l . QTI . Hence, at transition instants semi-Markov process behaves just 
as a Markov process, is used to combine the network topology (including radio channel) that 
are captured by the stochastic geometry with the routing protocol dynamics that are captured by 
the PGF framework. With semi-Markov process, the elapsed time between transitions may take 
several unit intervals, while at transition instants it behaves as a typical Markov process lf30ll . 
The semi-Markov process depends on the transition that occurred, though the transition instants 
follow a distinct probabilistic mechanism ||30|. Throughout our investigation we consider the 
states presented in Fig. [3] si, s 2 and s 3 corresponds to the transmission, relay and re-transmission 
states, respectively. The corresponding state transition matrix is given by 


P 


PSD (1 - Psd)psr (1 - Psd) (1 - Psr) 
1 0 0 


( 10 ) 


Psd (1 — Psd)Psr (1 — Psd){ 1 — Psr) 

where each element p tJ yields the transition probability from state i to j, as well as p SD and 
Psr yield the probability that the source message is correctly received by the final destination 
and relay node, respectively. 

Remark 3: The Markov process model has the property that a transition is made at every 
time instant. Differently, the semi Markov model is a more general class of processes where 
the elapsed time between transitions may take several unit time intervals. Then, it brings much 
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Fig. 3. Illustration of a Semi-Markov process with three states s;, i £ {1,2,3}. s i represents transmission, S 2 relay and S 3 
retransmission. 


more flexibility to the problem of modeling dynamic probabilistic systems. In fact, the successive 
state occupancies are governed by the transitions probabilities of an embedded Markov process, 
while the holding time in each state follows an integer-valued random distribution - the conflict 
resolution interval, see © in Section [IV] — which depends on both the current and the next states 


The PGF of the CRI Gn{z) characterizes the relay selection overhead at each hop. The 
probability that a certain CRI takes k slots can be retrieved from the respective PGF using the 
following proposition. 

Proposition 3: The probability mass function of L N , the length of the contention resolution 
interval (in transmission slots) conditional on the number of contending nodes is Pr[X — x] — 

0)/xl 


Proof: From the PGF definition it is straightforward that 


G n (z) = E[z x ] = ^}z x . (11) 

x=0 

Then, taking the nth derivative of (fTTI) and making z = 0 we attain the result in Proposition [3] ■ 
The holding time matrix incorporates into the embedded Markov chain of Fig. [3] the time 
spent in a given state between successive transitions. 



1 

1 + Ln 

1 + Ltv 

H = 

1 

1 

1 


1 

1 + Ln 

1 + Ln 


(12) 


where each element h t j yields the time elapsed in state i before transition to j. 
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A. Steady - State Throughput Efficiency 


Let i be an arbitrary state chosen as reference, and n, be its steady-state probability where the 
vector 7T corresponds to the stead-state distribution of the embedded Markov chain (see Fig. |3]). 
The time delay associated to the transition from state i to j is given by Dij. Similarly, R rj is the 
reward associated with the transition from state i to j. A cycle is defined as the time between two 
consecutive passages for the reference state states. A message is successfully delivered to the 
destination node when the process returns to the state ,s'i. The reward function R(r) corresponds 
to the cumulative reward after a period r and yields the total amount of correct receptions. 
According to the fundamental renewal-reward theorem lf30l . the long-term throughput is given 
by the following value, It is worth noticing that the resulting steady-state throughput efficiency is 
independent of the choice of the reference state. Moreover, the long term evaluation of R(r)/r 
actually corresponds to the normalized throughput of the protocol as introduced in ll32ll and 
revisited in If28l 


V 



T—> OO T 


E f=i fljRi 

EiiTTA’ 


(13) 


where K yields the number of states, R = p tJ R (/ is the mean rewarding, while H = 
pij h ij is the mean waiting time. 


VI. Numerical Results 

For both desired link and interferes, a composite fading channel model with Nakagami- 
m parameter of m = 16 (which corresponds to a Ricean fading channel with factor K = 
14.8dB) and shadowing standard deviation of 10 dB are considered as well as channel reciprocity 
between a communication pairs. Interferes are scattered over the network area with R rn = 25 
and Rm = 500 meters, while A = 5 x 10 _5 nodes/m 2 is considered (which corresponds to nearly 
40 interferers in average) and the path loss attenuation is set to 3. The CRI is conditioned on 
the mean number of candidate relays on the forward region of the source node. Without loss of 
generality, we consider three contending relays within source’s forwarding region to compute the 
relay selection interval. We assume that data and control information occupies one transmission 
interval. 

In a HD configuration, Fig. |4] shows the steady-state throughput efficiency fj for increasing 
separation distance between source-destination pair. As evidenced by the reactive relay curve 
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Fig. 4. Steady-state throughput efficiency for increasing separation distance between source - destination pair in HD 
configuration. 


(red dashed line), the steady-state throughput is severely compromised by the relay selection. On 
the other hand, when a fixed relay is considered within source’s transmission range, fj becomes 
much less susceptible to the degrading effect of longer separation distance between source and 
destination. 

Fig. |5] presents the steady-state throughput for increasing transmission power of source node, 
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Fig. 5. Steady-state throughput efficiency increasing source’s transmission power in HD configuration. 

whereas interferers transmit at 30 dBm. The separation distance between source and destination 
is kept at 50 m and the relay is randomly located in the source’s forwarding region. As expected 
the steady-state throughput improves with higher source’s transmission power. The reactive relay 
benefits most from high transmission power mainly because source can reach the destination more 
often without undergoing long relay selection intervals to select a suitable relay. 

The steady-state throughput efficiency for decreasing SI attenuation with FD configuration is 
shown in Fig.[6l The separation distance between source-destination is set to 50 m. As previously 
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Fig. 6. Steady-state throughput efficiency for decreasing SI attenuation in FD configuration. 

observed in the HD configuration, the steady-state throughput with reactive relay is severely 
degraded by the relay selection procedure and barely outperforms the HD scenario with 100 dB 
SI attenuation. It is also shown how the steady-state throughput efficiency scales for m = 2 and 
16. As can be seen, fj slightly degrades when the Nakagami-m parameter of the self-interference 
component increases from 2 (non line of sight) 16 (line of sight). 

Regarding the fixed relay scheme, the FD configuration with SI attenuation raging from 100 to 
nearly 70 dB show much better performance than HD mode, whereas the performance degrades 
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faster with lower attenuation values. From 120 to 100 dB the FD gain is negligible, because SI is 
low and the source node reaches the destination directly. However, for lower SI attenuation the 
source node uses the relay more often which increases fj (the reward doubles). Despite that, if 
the attenuation is lower than 70 dB, the SI also compromises the relay performance and degrades 
the Steady-state throughput efficiency considerably. Moreover, as indicated by the dot-dashed 
threshold lines, the FD configuration outperforms the HD one for high SI attenuation values. 
However, if the SI attenuation drops below 70 dB, the FD performance degrades and HD becomes 
a better alternative. In fact, it would be advantageous to consider a dynamic forwarding scheme 
by which potential relays would change between FD and HD depending on the attainable SI 
value. 


VII. Conclusions and Final Remarks 

In this work, we assess how in-band full duplex relaying performs for distinct selection 
strategies, namely fixed and reactive. Our investigations are carried out using an analytical 
framework based on stochastic geometry and semi-Markov process. The former is used to model 
the network deployment area and radio channel, while the latter characterizes the cost of the 
relay selection procedure under study. Our results show that the benefits of the FD configuration 
are compromised by the frequent selection of the next hop relay. Moreover, the steady-state 
throughput degrades depending on the relative distance between source, relay and destination, 
as well as the respective SI attenuation. As pointed out in lfl4l . Ifl5l . a binary search tree 
implementation does not take full advantage of the side information (local network topology) 
to shorten the relay selection time. Hence, we plan to extend this work by considering Dutch 
auction-based relay selection to exploit the local topology of the network in order to shorten 
the selection transactions. We also expect to assess the secrecy of the network, as in lf33l . by 
employing a similar framework. 
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